Entanglement Evolution in a Non-Markovian Environment 
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We extend recent theoretical studies of entanglement dynamics in the presence of environmental 
noise, following the long-time interest of Krzysztof Wodkiewicz in the effects of stochastic models of 
noise on quantum optical coherences. We investigate the quantum entanglement dynamics of two 
spins in the presence of classical Ornstein-Uhlenbeck noise, obtaining exact solutions for evolution 
dynamics. We consider how entanglement can be affected by non-Markovian noise, and discuss 
several limiting cases. 
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I. INTRODUCTION 

A quantum system of interest may be large or small, 
but its background, called an environment, is almost 
always complex, and is often represented by a bath of 
bosons or fermions, or by classical random fields. In all 
these cases, system dynamics is described by a quantum 
master equation that governs the evolution of the reduced 
density matrix of the system. In the current decade, de- 
coherence dynamics of entangled quantum systems un- 
der the influence of environmental noises has been ex- 
tensively discussed in different contexts involving atoms, 
ions, photons, quantum dots, and Josephson junctions, 
to name several. This is all related to new regimes of 
information processing, such as quantum cryptography 
and quantum computation [l[ . An important category of 
such research has treated the fascinating domain where 
entanglement of qubits evolves even though the qubits do 
not interact, even indirectly. An example is sketched in 
Fig. [T]and we restrict our attention here to this category. 

In experimental environments an entangled system 
may be exposed to vacuum noise, phase noise, thermal 
noise, and various classical noises, as well as mixed com- 
binations of noises. A number of idealized models have 
provided new insights by allowing entanglement evolu- 
tion to be followed by solving the appropriate quantum 
master equationfsee Zyczkowski, et al. Q, Daffer, et al. 
0| , as well as 0, H, H 0,11, [§] ) • Most research on entangle- 
ment dynamics has been focused on ambient noises from 
environments that obey the Markov (no memory) as- 
sumption. Recently, there is growing interest in the non- 
Markovian entanglement dynamics for both discrete and 
continuous quantum systems (see [H EH EH 0, EH EH 



and the overview in Q). 



In truth every environment is non-Markovian. Non- 
Markovian noise was a repeated theme in the research of 
Krzysztof Wodkiewicz [ij, EH El, El M, SH H pi pH 
HI HI H3| > and we present our findings as a contribution 
to his scientific memory. 

As far as we know there are no fully systematic in- 
vestigations of non-Markovian noises or of their effect 
on the coherence dynamics of non-interacting spin sys- 
tems. In particular, a perturbative theory leading to the 
Markov approximation is still lacking. The purpose of 
this paper is to present a study of such problems in the 
simplest form. We will consider classical non-Markovian 
noises, modelling them as so-called Ornstein-Uhlenbeck 
processes, and derive the consequences for entanglement 
dynamics. This can be considered an extension of our 
earlier note on entanglement sudden death (ESD) under 
classical Markov noises [28| . 
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FIG. 1: Sketch of remote qubits in a quantum memory net, 
'Email address: ting.yu@stevens.edu where the dashed lines indicate entanglement, but not inter- 

tEmail address: eberly@pas.rochester.edu action. 
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II. THE KUBO-ANDERSON MODEL 
EXTENDED TO TWO QUBITS 

We consider an entangled pair of spins both of which 
are subject to frequency fluctuations that are random 
[29l | . We adopt a model for these fluctuations that 
treats them as caused by noisy environments described by 
Ornstein-Uhlenbeck processes. This well-known Gaus- 
sian noise model is non-Markovian in the general case but 
has a well-defined Markov limit. To focus exclusively on 
the effects on the entanglement of the spins as it arises 
from the noise, we assume the spins to be affected sepa- 
rately by separate environments, and not to interact with 
each other in any way, especially not through the noises. 
Thus the spins could be, for example, remote compo- 
nents of a quantum memory net (as in Fig. [T|) under 
steady attack by weak local noise. This compromises the 
preservation of their entanglement. 

The Hamiltonian of the two-spin system can then be 
formally written as (we set ft = 1 ) : 



«toH*J - — ^ — a z H ^ — < 



(1) 



where fi>i(t) and f2_e(£) are the independent fluctuations 
of the spin transition frequencies (level spacings). They 
have the mean value properties 



M[Qi{t)} 

M[su(t)su{s)] 



0, 

a(t — s) 



r n _ 



j\t-s\ 



= A,B, 



(2) 



(3) 



where M [•] stands for the statistical mean over the noises 
£lA(t) and fis(t). Note that 7 is the noise bandwidth, 
and 7 _1 = r c defines the environment's finite correlation 
time of the noise. For simplicity, we will take the noise 
properties to be the same for A and B (e.g., Ta = = 
r), although independent. In the limit 7 — ► 00, Ornstein- 
Uhlenbeck noise reduces to the well-known Markov case: 



a(t, s) = TS(t — s). 



(4) 



For the total system described by the Hamiltonian (JTJ) , 
the stochastic Schrodinger equation is given by 



ij t m)) = H to t(t)\m)- 



(5) 



The explicit solution for the stochastic Schrodinger equa- 
tion can be readily obtained in terms of a stochastic uni- 
tary operator: 

m)) = u(t,n A ,n B )\i;(o)), (6) 

where the stochastic propagator U(t, Ha, ^b) is given by 

U(t,Q A ,n B ) = e -^o d < nA ^ +UB ^). (7) 

The reduced density matrix for spins A and B is then 
obtained from the statistical mean 



P (t) = M[\m)m)\] 



(8) 



The master equation for the reduced density matrix for 
the two-spin system in a non-Markovian regime can be 
readily derived from the stochastic Schrodinger equation 



dp 
~d~b 



G(t) 



(2p 



A A 
°z P°z 



B B\ 
°z P°z ) 



where 



G{t) = j a(t-s)ds=^(l 



,-7* 



), 



(9) 



(10) 



The memory information of the environmental noises is 
encoded in the time-dependent coefficient G(t) where 
t c = I/7 characterizes the environmental memory time. 
In the Markov limit r c -> (7 -> 00), when G{t) -> T/2, 
equation © reduces to the well-known Markov master 
equation in the presence of dephasing noises. 



III. EXACT SOLUTIONS FOR QUANTUM 
EVOLUTION 



Solutions of master equations for the noisy evolution of 
two-spin density matrices in terms of the Kraus operator- 
sum-representation have been given before (see, for ex- 
ample, [H, [Hi H3]). In many cases of physical interest, 
the Kraus representation allows a transparent analysis of 
entanglement decoherence without invoking the explicit 
forms of the initial conditions. In what follows, we will 
use the fact that for any two-spin initial state p (pure or 
mixed), the evolution of the reduced density matrix can 
be written compactly as 



p(t) = J2K^t)p(0)Kl(t), 



(11) 



where the Kraus operators satisfy KjJCft = 1 for 
all 

In order to derive the desired Kraus operators for the 
reduced density matrix we begin by noting that the so- 
lution for just spin A can be written: 



where 



\m) = u(n A ,t)mo)) 



U{fl A ,t)=exp [-iF(t)<r z ] 



(12) 



(13) 



with the stochastic process F(t) = f*dsn A (s). Then 
our first task is to express the stochastic density operator 
Pst = I'^WXV'WI m the Kraus- like operator representa- 
tion form: 



p st {t) = exp [-iF{t)cr z ] p(0) exp [iF(t)a z 



(14) 



where p(0) — \t/j(0)) (ip(0)\ is the initial state of the sys- 
tem, which is assumed to be independent of the noise. 
The desired Kraus operators for the spin are obtained by 
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taking a statistical mean over the noise f2^i(f) for qubit 
A and are given by 



Et = 



PA{t) 



, E 2 = 



q A {t) 




1 

where the time-dependent Kraus matrix elements are 



(15) 



QA(t) 
P A {t) 

/(*) 



l-p\(t), and 
exp [-/(*)], with 

G{s)ds 
2 7 



(16) 
(17) 



(18) 



and similar expressions for ps(t) and <?B(f)- The two- 
qubit case given here can be easily applied to N noninter- 
acting qubits, an extension we reserve for later attention. 

Since our two spins are evolving independently, we 
have the following four Kraus operators in terms of the 
tensor products of E\ and E 2 : 



K x = 
K 2 = 
K 3 = 
K 4 = 



PA 



PA 



q A 


q A 




PB 

1 

q B 


PB 

1 

q B 




(19) 
(20) 
(21) 
(22) 



IV. 



NON-MARKOVIAN ENTANGLEMENT 
DYNAMICS 

A. X matrix and concurrence 



We now consider entanglement dynamics of two half- 
integral spins (qubits) with an initial density matrix with 
the common X-form 1311: 



„ab 



(pix pi4 N 

P22 P23 

P32 P33 

V/941 ,044/ 



(23) 



Such X states occur in many contexts and include pure 
Bell states as well as Werner mixed states. 

For two qubits, entanglement can be evaluated un- 
ambiguously via the concurrence function [33| . which 
may be calculated explicitly from the density matrix 
p AB . For qubits A and B we have: C AB = C{p AB ) = 
max{0, Q(t)}. Here Q(t) is defined as 



Q 



A 4 , 



(24) 



where the quantities Xi are the (generally time- 
dependent) eigenvalues in decreasing order of the follow- 
ing (nonlinear in p) matrix: 



C = pK^<V>>^), 



(25) 



where p* denotes the complex conjugation of p in the 
standard basis +, +), |+, — ), — , +), | — , — ), and a y is the 
usual Pauli matrix expressed in the same basis. 

From the general solution p9H22[) . one can easily show 
for the initial state (l23l) that one finds 



Q(t) = 2max^|p 3 2(i)| - vW°W(0), 

|Pl4(*)| - V P22 (0)/O 33 (0) 



(26) 



B. Solutions for non-Mar kovian disentanglement 

The Ornstein-Uhlenbeck phase-noise solutions for the 
density matrix elements of a general initial state are given 
by 



P12(*) 


- Pi 2 (0)e-W 


(27) 


Pl3(t) 


= Pi 3 (0)e-^W, 


(28) 


m(t) 


= pu(0)e~ m , 


(29) 


Pu(t) 


= p 3i (0)e~ f{t) 


(30) 


P23(*) 


= P2s(0)e-^\ 


(31) 


Pu(t) 


= p u {0)e- m) 


(32) 


Pii(t) 


= Pii (0) (i= 1,2,3,4) 


(33) 



where f(t) is defined in (JTSJ). Let us note that in the limit 
7 — > oo, we recover the standard Markov approximation 
where f(t) = Tt/2. 

Although there is no compact analytical expression for 
the concurrence C{p(t)) with an arbitrary initial state, 
we can readily show that preservation of entanglement is 
restricted by the inequality 



C(p(t)) < e-WC(p(0)). 



(34) 



A sharper result occurs for the X matrix under consid- 
eration because the diagonal elements are independent of 
t. Since they all vanish as exp(— /(f)) for increasing f , we 
know that Q{t) must eventually become strictly negative 
if diagonal values are initially non-zero (e.g., any finite- 
temperature equilibrium state). Negative Q mandates 
C AB = 0, so ESD must occur. Next we will consider key 
limiting cases. 



C. Entanglement decay: Stationary limit 

We consider now the stationary limit 7* ^> 1. Then 

r 



/W-&< + -(^-l)] 
2 7 



-f 



(35) 
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FIG. 2: The graph shows C AB vs. Ft and 7/r. The value 
a = 1/3 has been chosen. The reservoir bandwidth 7 controls 
the approach to the Markov limit, and we see that for 7 < V 
the inevitable onset of ESD (the region where C(t) — 0) can 
be substantially delayed. 



Therefore, from we get 

C( P {t)) < e- rt C(p(0)). 



(36) 



Hence, the entanglement decay rate is at least as rapid as 
r, and may be much faster. Clearly, the stationary limit 
is identical to the Markov limit a(t — s) = TS(t — s). 



D. Entanglement decay: Short-time limit 

Now let us turn to the opposite limiting case: jt -C 1. 
In this case, we can use the following approximation, 



Therefore, 



^~l- 7 t + i 7 2 ; 2 



f(t) = i r 7 t 2 



Then concurrence decay is bounded by 
C(p(t)) < e^ t2 C(p(0)) 



(37) 



(38) 



(39) 



Hence, the effective disentanglement time is given by: 



Tdii 



(40) 



Clearly, for non-Mar kovian noises, the short-time limit 
is more interesting since it shows that the resultant en- 
tanglement behavior deviates significantly from the well- 
known Markov dynamics. Obviously, the smaller 7 is, the 
better approximation we have. For both limiting cases, 
it is easy to prove that a sufficient condition for ESD to 

OCCUr is ,011^22/333/044 ^ [H]. 



V. EVALUATION FOR A SPECIAL X STATE 

In Fig. [2] we show the evolution of the concurrence for 
a specific X-form entangled state: 



P a (0) 



fa 
Oil 
Oil 

\0 l-aj 



(41) 



where < a < 1, so the initial concurrence is C(0) = 
2/3[l - y/a(l - a)} > 0. 

The time dependence of the concurrence of this state 
is well-known in the Markov dephasing limit [28| . For 
our present non-Markovian case, which introduces the 
environmental bandwidth 7, the time dependence of the 
Q parameter satisfies 



Q(t) = ~(e-/«- A /^T^)). 



(42) 



For all finite values of 7 the quantity e~^ l > approaches 
zero exponentially at long times, and then Q(t) must 
become negative, so ESD inevitably occurs, with a finite 
disentanglement time t^su given by 



,-2/(t ESD ) _ 



a(l — a). 



(43) 



Fig. [3] shows several interesting features of entangle- 
ment evolution under Ornstein-Uhlenbeck noise. Clearly, 
we see that non-Markovian noises have markedly different 
affects on entanglement evolution at short times, while 
the long time limit gives rise to familiar Markov behav- 
ior. First, note that ESD must occur for the entire pa- 
rameter range of a, except for the end points a = 0,1, 
given our initial X-state. However, the ESD times are dif- 
ferent for non-Markovian short-time and stationary limit 
cases. From Fig. 3a, it should be noted that the dis- 
entanglement time for non-Markovian regimes could be 
significantly longer than the disentanglement times in the 
Markov limit if the dissipation is small. However, once 
the state becomes separable it will never become entan- 
gled again. That is, entanglement rebirth or revival does 
not occur for Ornstein-Uhlenbeck noise UM. 



VI. CONCLUDING REMARKS 

We have presented here, as a contribution to the sci- 
entific memory of Krzystof Wodkiewicz, the first results 
of a new investigation that has clear connections to his 
long-time interest in quantum systems evolving under the 
influence of stochastic perturbations. One of the targets 
of his creativity and energy, over many years, was the 
challenge presented by the influence of non-Markovian 
noise and we have addressed that challenge with some 
calculations focused on entanglement. 

We have shown that entanglement dynamics under 
Ornstein-Uhlenbeck noise can be affected in several dif- 
ferent ways, depending on the initial entangled states and 
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the noise correlation time. It can be seen that the non- 
Markovian properties can prolong the life of entangle- 
ment. We note that the effective long-time relaxation 
rate T is ordinarily associated with experimentally ac- 
cessible relaxation times such as T\ and T^. Fig. [2] high- 
lights the unusual domain 7 < T, in which these times 
are shorter than the internal environmental relaxation 
time t c . Entanglement survival is of fundamental inter- 
est at short times in quantum information processing (see 
dHHl]). In the case of the short-time limit our results 
capture the features of quadratic rather than exponen- 
tial decay at early times. In this simple model, non- 
Markovian noises appear to play a role as a short-time 
decoherence buffer, but entanglement measured by con- 
currence will inevitably conform to the stationary limit 
at long times. Finally, it is interesting to note that our 
findings based on the classical phase noise model can be 
extended into the case of quantum phase noises where the 
environment is modeled as a set of harmonic oscillators 
at a finite temperature [37j . 



FIG. 3: The graphs show C^ B vs. t with a = 1/3. Fig. 
3a shows that with j/T = 0.5, an initially non-Markovian 
entanglement (solid line) evolves in a markedly different way 
compared to its Markov limit (dotted line) . Clearly, the short- 
time limit (dashed line) gives a better approximation than the 
stationary limit (dotted line). However, as shown in Fig. 3b 
where 7/r = 5, the difference is washed out at later times 
when the short-time limit (dashed line) ceases to be a good 
approximation. 
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